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ABSTRACT 

We model the effects of sliocks on the diffuse, X-ray emitting baryons in 
clusters of galaxies. Shocks separate the infalling from the inner gas nearly at 
equilibrium, and dominate the compression and the density gradients of the 
latter in the dark- matter potential of the cluster. We find that, independently 
of the detailed shape of the potential, the density gradient is steeper and the 
compression factor larger for the richer clusters. We show, considering the 
different merging histories, that in the hierarchical cosmogony the above effects 
lead, in X-rays, to a luminosity-temperature relation L oc at the scale of 
groups which flattens down to L oc for rich clusters in accord with the 
observations, and then saturates toward L oc for higher temperatures. Prom 
the merging histories we also compute statistical fluctuations of the L — T 
correlation. 

Subject headings: galaxies: clusters - galaxies: intergalactic medium - 
hydrodynamics 
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1. Introduction 

The X-ray emission from clusters of galaxies enables direct probing of gravitationally 
bound and virialized regions with virial radii of a few Mpcs, comprising total masses 
M ~ 10^^ Mq mostly in dark matter (DM hereafter). 

On the one hand, the X-ray temperature T oc GM/R^ measures the depth of the 
potential wells. On the other hand, the bolometric luminosity L (x n R\ T^l'^ emitted as 
thermal bremsstrahlung by the intracluster plasma measures the baryon number density n 
within the volume The L — T relation constitutes a crucial link between the physics of 
the baryon component and the dynamical properties of the DM condensations. 

The simplest model describing the former holds n to be proportional to the average 
DM density p, and Rx to i?^, so that n oc p oc M/Rf, obtains (self-similar model, hereafter 
SS, Kaiser 1986). If so, the luminosity would scale as L oc p^/^T^, which is inconsistent 
with the observed correlation close to L oc (Edge & Stewart 1991; Mushotzky 1994; 
Tsuru et al. 1996). Further steepening at the temperatures of galaxy groups is indicated 
for the emission not associated with single galaxies (Ponman et al. 1996). In addition, 
the SS model when combined with the standard hierarchical cosmogony (see Peebles 1993) 
yields for the clusters a local X-ray luminosity function too steep or too high compared 
with the data (Evrard & Henry 1991; Oukbir, Bartlett & Blanchard 1996). 

So the indication is that the ratio n/p is to depend on M or T. It will be 
convenient to write the volume-averaged in terms of two factors: the compression 
factor g{T) > 1, describing the gas overdensity relative to the outer value at the 
"boundary", taken here to be at Ry as discussed in §3; and the inner shape factor 
I{Ry,T) = R~^ (fr n'^{r,T)/n^{Ry,T), with the main contribution coming from inside 
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Rx- The result is: 

L oc g\T)I{R,,T)RlT'/'p' oc g\T) I{R,,T) p^/^ , (1) 

where the last term follows from expressing R^ from T oc M/R^ ex pRl- To obtain the 
average L — T relation to be compared with data, the factor g'^{T) in eq. (1) has to be 
averaged over the cluster histories, as we carry out in §2.3. 

So, the difference of eq. (1) from the SS model has been factored out into the terms 
g'^{T) and I{Ru,T), which are determined by the hydro- and thermodynamics of the gas in 
the forming cluster wells. At ^ ;^ 1 the gas is expected to be preheated by feedback effects of 
star formation: injections of energy of stellar origin like Supernova winds eject the gas from 
the shallower potential wells, and heat the residual and the ejected gas to temperatures 
Ti ^ 10^ K (Dekel & Silk 1986; Kaiser 1991; Ciotti et al. 1991; David et al. 1993, 1995; 
Cavaliere, Colafrancesco & Menci 1993). In addition, preheating is necessary to prevent 
too short cooling times (as pointed out by Cole 1991, and addressed by Blanchard, Vails 
Gabaud & Mamon 1992; see also White & Rees 1978) in early potential wells, shallower on 
average but containing denser gas. Subsequent evolution will lead to an increasing recovery 
of the universal baryonic fraction (White et al. 1993). 

Previous attempts (Kaiser 1991; Evrard & Henry 1991) to tackle the L — T relation 
are based on the extreme assumption that the gas inside the X-ray core Rx is preheated but 
never subsequently shocked or mixed. Here we discuss the o^/ier extreme, i.e., the effects of 
shocks and mixing on the gas density inside clusters. 

2. The Shock Model 

As in the collapses the gas velocity becomes supersonic, shock fronts form at about 
Ry, and separate the infalling from the inner gas already at virial temperatures. In fact. 
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numerical simulations (see Evrard 1990; Takizawa & Mineshige 1997) of isotropic collapses 
show that, when the outer gas temperature is appreciably lower than the virial value T, a 
spherical shock front forms and, in the vicinity of Ry, slowly expands outwards leaving the 
gas nearly at rest and with a nearly fiat temperature profile. When realistic, anisotropic 
collapses are considered (Navarro et al. 1996, Tormen 1996) shock fronts still form and 
convert into heat most of the hydrodynamical energy (Schindler & Miiller 1993; Schindler 
& Bohringer 1993; Rottiger, Burnes, & Loken 1993). 

Across the shock the gas entropy rises, and correspondingly a jump in the gas density 
and in the temperature is established from the exterior values rii, Ti to the interior ones 
n2, T2. The density jump provides a boundary condition for the inner gas distribution in 
the form of the compression factor g{T2/Ti) = n2/ni. In addition, the interior temperature 
T2 governs, at equilibrium in a given gravitational potential, the inner density profile and 
hence the shape factor I{Rv, T). These two effects enter eq. (1) for L, and will be discussed 
in turn. 



2.1. The Compression Factor 

The values of n2, T2 and of the interior gas velocity V2 are related to their outer 
counterparts ni, Ti and Vi by the requirements of mass, momentum and energy conservation 
across the shock. The plasma behaves as a perfect gas with three degrees of freedom, and 
the Hugoniot adiabat (see Landau & Lifshitz 1959) yields the compression factor 

1/2 



For strong shocks with T2 » Ti, this saturates to the value g — ^, while for T2 — > Ti it 
attains its lowest value g — 1- 

The pre-shock temperature 71 is provided by the stellar (thermonuclear) energy 



-6- 



feedbacks recalled in §1, or by the virial (gravitational) temperature inside the clumps 
which are to merge with the cluster during its merging history considered below. In fact, 
the stellar feedbacks set for Ti the lower bound Tu which we identify with the lowest 
temperatures (around 0.5 keV, Ponman et al. 1996) measured in groups. 

Pre-shock temperatures of this order do not affect the rich clusters; instead, they 
affect the compression factor in the shallower potential wells with T ~ 1 keV, as prevail at 
redshifts z ^1 but are also present at z ~ 0. The full behavior of g{T2/Ti) in shown by the 
dashed line in fig. 1 when Ti = Ti^, and the latter is in the range 0.5 — 0.8 keV. The actual 
values of Ti will be discussed in §2.3, taking into account the merging histories. 



The post-shock temperature T2 can be calculated from the pre-shock velocity Vi. The 
latter is driven by the gravitational potential V{r), and reads Vi — [—aV{Rv)/mp]^^'^ with 
a = 2[1 — V{Rm)/V{Ry)]. Here Rm is the radius where infall becomes nearly free; its upper 
bound is obtained by equating the Hubble flow to the free-fall velocity, which yields a ~ 1.4. 
The post-shock condition is closely hydrostatic, i.e., V2 <^ vi, as shown by the simulations. 
Then, using the equations in Landau & Lifshitz (1959) with the above value of fx, we find 



The inner gas profile relative to that of DM, p(r) say, is governed at equilibrium by the 
scale-height ratio /5 = iimpa^/kT2, where /i ~ 0.6 is the gas mean molecular weight, m.p the 
proton mass, and (7^ the one-dimensional velocity dispersion of the DM. The profile 



2.2. 



The Gas Disposition 



kT2C^~V{R,) + lm . 

o o 



(3) 



n{r) oc [p(r)]^(^) 



(4) 



(Cavaliere & Fusco-Femiano 1976), applies for a nearly fiat T{r) const T2. 



- 7- 



The function P{T) entering the profile (4) is easily computed from eq. (3), for a 
given DM potential V(r) corresponding to p(r). For the King potential (see Sarazin 1988) 
V{r) = —9iJ,mpa^rcln[{r/rc) + (1 + r^/r^)^/^]/r with the core radius Tc = Rv/12, we obtain 
P{T) increasing somewhat from the value P 0.5 for T Ti to /3 fa 0.9 for T > Ti. A 
similar result obtains using the potential proposed by Navarro et al. (1996). These two 
instances are illustrated in fig. 2, and demonstrate that in all cases the static gas density 
profile is shallower for lower T clusters. 

2.3. L — T from Merging Histories 

For a given Ti the compression factor g is computed after eq. (2), and the shape 
factor I{R,T) is obtained by integration of n^{r) computed after eq. (4); then the L — T 
relation may be obtained from eq. (1). But Ti depends on the thermal conditions of the 
infalling gas, which is preheated by Supernovae and further heated through virialization 
inside merging clumps. 

In the former case, we take Ti = Tu — 0.5 — 0.8 keV with a flat distribution, and 
obtain for the dashed line in fig. 1. In the latter case, repeated merging events introduce 
fluctuations of Ti above Ti*, and hence of the interior density n2, which modify g^. The 
average effect is shown by the solid line in fig. 1, and the corresponding variance is 
illustrated by the shaded area. 

To include both conditions, we take Ti to be the higher between the preheating value 
Ti* and the virial temperatures prevailing in the clumps accreted by the cluster. We perform 
a statistical convolution of the L — T relation over such merging histories, using Monte 
Carlo realizations of hierarchical merging trees of dark halos as introduced by Cole (1991). 
The code, written by one of us (P.T.), is based directly on the excursion set approach of 
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Bond et al. (1991) to the mass distribution. 

For each merging event concerning a cluster with a current virial temperature T', 
we compute g{T'/Ti) for the clumps being accreted, weighted with the associated mass 
fraction. We show in fig. 1 the quantity {g'^) averaged over the histories ending into a 
cluster of temperature T, along with its dispersion. These two quantities are used in eq. 
(1) to predict the average L — T correlation and its scatter. The results are shown in fig. 
3, for the DM potential of Navarro et al. 1996; a different V{r), like King's, only steepens 
somewhat the low-T behavior. 

3. Results and Discussion 

Here we have proposed a model for the intracluster gas, to capture in a simple way 
one essential component of the complex gravitational systems constituted by groups and 
clusters. The model is focused on the formation of shocks between the the gas nearly in 
equilibrium with the cluster potential, and the infalling one. The latter is preheated by 
the release of thermonuclear energy, or by the gravitational energy in subclusters. Shock 
heating and compression determine the L — T relation. 

The latter comprises both clusters and small groups in a single dependence, which 
smoothly flattens from L oc (for groups with T ^ 2 keV), to L oc T^-^ (for clusters 
with 2 keV^ T ^7 keV), toward L (x T'^ (for higher T). Such behavior fits both the 
cluster data (Edge & Stewart 1991) and those for groups (Ponman et al. 1996) which - 
if considered separately - would require a much steeper L — T relation. Correspondingly, 
the volume-averaged baryonic fraction grows by a factor around 3 from small groups to rich 
clusters, but remains within 1.3 times the universal value. 

In addition, we predict the gas density profiles to have a fiat central region (the gas 
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"core", see fig. 2), independently of the detailed shape of the DM potential; the profiles are 
actually steeper for larger virial temperatures. Correspondingly, the size Rx of the X-ray 
core (defined, e.g., at one half the integrated emission) grows with mass slower than M^/^. 
The average cooling time within Rx exceedes the Hubble time out to 2; fa 2, differently 
from the SS model. 

We have checked that these results persist when we relax the approximation of a flat 
temperature profile in the cluster, and adopt instead a polytropic distribution with index 7 
ranging from 1 to 5/3; for For 7 > 1 the temperature declines from the center toward the 
shock position. 

In time, the shocked region expands and outgrows the infall velocity decreases, 
and the shock weakens with T2 approaching Ti. However, this occurs only over several 
dynamical times as shown by the N-body simulations (see Takizawa & Mineshige 1997); 
meanwhile, the shock positions remain close to as taken here. 

We do not stress, instead, the 2;-dependence (1 + z)^'^ of the normahzation provided by 
the factor p^^'^(z) appearing in cq. (1). In fact, such dependence is easily swamped by the 
place-to-place variations of rii and by the systematic increase in contrast of the large scale 
structures hosting groups and clusters (see Ramella, Geller, & Huchra 1992). 

The robust predictions of the shock model do not require spherical symmetry, but 
only a small bulk velocity of the inner gas compared to {GM/rUp R^Y^^. Thus the model 
includes anisotropic, recurrent merging with other clumps of dark and baryonic matter. 

The effects of extreme merging events are as follows. The few events involving 
comparable subclusters reshuffle the baryonic content and mix its entropy, but only 

moderatly affect temperature and density. At the other extreme, the more isotropic 
accretion to a cluster of many small condensations gravitationally heated at temperatures 
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Ti 1 keV yields the highest compression and the largest contribution to the X-ray 
luminosity. Our Monte Carlo simulations span the range between these extremes. 

The stellar preheating at Ti* of the external gas is essential to provide a lower hmit for 
Ti. This ensures that the accreting gas, even when in a very shallow well or in a diffuse 
state, starts on a relatively high adiabat, corresponding to Ti^. = 0.5 — 0.8 keV. Such a 
heating of thermonuclear origin breaks down the otherwise self-similar form n/p =cost of 
the ratio of gas to DM density to yield, at the boundary, the form n/p (x g{T/Ti) shown in 
fig. 1. Values of Ti^, smaller than 0.5 keV would lead, in a strictly self-similar evolution of 
DM halos, to L (X at variance with the data; larger values to a severe depletion of the 
gas and of the luminosities in groups and clusters. 

The merging histories also produce the considerable scatter in the L — T relation 
shown in fig. 3, since the different virial temperatures of the stochastically merging clumps 
induce intrinsic variance in the internal density n2. Further scatter may be contributed by 
the vagaries in the ambient density rii , and by the possible lack of dynamical equilibrium 
in some groups, see discussion by Governato, Tozzi & Cavaliere 1996. 

The shock model in the simple form presented here applies to the gas settled to 
equilibrium after each dynamical perturbation. This takes sound propagation times, 
somewhat shorter than the dynamical timescale taken anyway by the DM to adjust to 
equilibrium (Tormen 1996). The residual converging motions, even in spherical N-body 
simulations (Takizawa & Mineshige 1997), tend to balance the expansion of the shocked 
region to yield only small net velocities V2 ~ 100 km/s. These may be associated with some 
adiabatic compression of the central regions, but the resulting heating is only mild, as long 
as shocks form at radii of order R^. 

The other extreme is tackled by the model proposed by Kaiser (1991) and refined by 
Evrard & Henry (1991). This assumes that, after preheating, the central cluster region 
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visible in X-rays contains the same gas (about 10 % of the present total) engaged in a 
smooth adiabatic compression. However, we find only 5% of the largest single progenitors 
to have masses (DM and hence gas) exceeding 10% of the present values at ^ > 1.5, when 
most stellar preheating takes place. On the other hand, N-body simulations (see refs. in §2) 
and observations (see Zabludoff & Zaritsky 1995) show that each cluster history includes a 
few merging events between comparable structures, which will reset the core gas to a higher 
adiabat. 

The adiabatic model (with preheating) predicts a single, scale-free relation L oc T^-^, 
OT L (X if the gas equilibrium holds out to Ry as in Evrard & Henry (1991). However, 
on the largest scales that ought to sample fairly the universal baryonic fraction (White et 
al. 1993) one expects saturation toward the scahng L oc of the SS model; in addition , 
at the group scales a much steeper dependence is indicated by observations. Such opposite 
departures of the L — T correlation from a single power-law are beyond the reach of the 
adiabatic model, but within the predictions of the shock model. 

The model we propose leads (Cavahere, Menci, Tozzi 1997) to a local luminosity 
function N{L, 2; = 0) in agreement with the data, and to N{L, z) with the mild or no 
evolution shown by recent data, and confirmed by the deep X-ray counts. 

We acknowledge informative discussions with M. Ramella, the helpful comments of the 
referee, and grants from MURST and ASI. 
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Fig. 1. — The square compression factor g^{T/T-i) is shown by the dashed hne when Ti 
takes on its lowest value Ti*, with the latter uniformly distributed in the range 0.5 — 0.8 keV. 
The different values of Ti due to the merging histories (discussed in §2.3) affect the average 
dependence as shown by the solid line, and provide the 2-(t dispersion shown by the shaded 
region. A tilted CDM power spectrum in a critical cosmology (see White et al. 1996) has 
been used. 

Fig. 2. — The gas density profile (for a uniform T) derived from eqs. (3) and (4) using the 
King DM potential (upper panel), or (lower panel) that given by Navarro et al. (1996). The 
dotted lines refer to a group with T — Ti —0.8 keV, and the solid lines to a rich cluster with 
Ti < T = 8 keV. 

Fig. 3. — The L-T relation from the shock model, convolved with the merging histories 
of DM halos, is compared with data for clusters of galaxies (filled squares, from Edge & 
Stewart 1991), and for groups (open squares, from Ponman et al. 1996). The DM potential 
of Navarro et al. (1996) is used. The shaded region outhnes the 2 — a scatter expected from 
the merging histories. 
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